Partially asymmetric exclusion models with quenched disorder 
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We consider the one-dimensional partially asymmetric exclusion process with random hopping 
rates, in which a fraction of particles (or sites) have a preferential jumping direction against the 
global drift. In this case the accumulated distance traveled by the particles, x, scales with the 
time, t, as a; ~ t^^^ , with a dynamical exponent z > 0. Using extreme value statistics and an 
asymptotically exact strong disorder renormalization group method we exactly calculate, Zpr, for 
particlewise (pt) disorder, which is argued to be related as, Zst = Zpr/2, for sitewise (st) disorder. 
In the symmetric case with zero mean drift the particle diffusion is ultra-slow, logarithmic in time. 

PACS numbers: Valid PACS appear here 

Driven lattice gas models are able to describe the 
generic features of stochastic transport far from equilib- 
rium Jj, ^, 3]. This kind of transport is often observed, 
e.g. in transport systems 3 or (micro-) biological appli- 
cations [sj , where a steady flow is maintained by a steady 
input of energy. As an important recent application we 
mention the bi-directional transport on microtubules , 
in which motion of cargos is realized by attaching a few 
molecular motors of different type. In our approach a 
cargo-motor complex is modeled as an effective parti- 
cle with random hopping rates, the value of which de- 
pends on the number of attached motors. This kind of 
molecular machine allows for an efficient navigation in 
the cell, because the velocity can be controlled via the 
typical number of attached motors. Moreover these mod- 
els can be related via simple mappings to other problems 
of nonequilibrium statistical physics which include e.g. 
surface growth problems jll Q ■ 

Here we discuss the most prominent model of this 
kind the so-called asymmetric simple exclusion process 
(ASEP), what is considered on a periodic chain of N sites 
and with M particles. For particlewise (pt) disorder par- 
ticle i may hop to empty neighboring sites with rates pi 
to the right and qi to the left, where pi and qi are in- 
dependent and identically distributed random variables. 
For sitewise (st) disorder the random hopping rates are 
assigned to given sites of the lattice. In the totally asym- 
metric model qi = 0, for every i, whereas in the partially 
asymmetric model, qi > pi for a finite fraction of i. 

A characteristic feature of the ASEP and its variants is 
their sensitivity to spatial inhomogeneities or quenched 
disorder of any kind. This is expressed in boundary 
induced phase transitions or phase separated states 
caused by a single defect or disordered lattices Q . For 



pt disorder some exact results are available 0, 0, ^2 , 
in particular for the totally asymmetric case. Depend- 
ing on the extremal properties of the hopping rate dis- 
tribution there is a dynamical phase transition in the 
system, separating a homogeneous state from a nonho- 
mogeneous one, in which there is a macroscopic parti- 
cle free region in front of the domain of occupied sites. 
For st disorder the analytical results are scarce 

PI 

and 

most of our knowledge is based on Monte Carlo simula- 
tions and mean- field calculations 9, 15]. In the 
totally asymmetric model as the strength of disorder is 
increasing in the stationary state the particle density is 
changing from a homogeneous phase into a segregated- 
density phase, in which macroscopic regions with den- 
sities 0<pc<l~Pc<l coexist and the macroscopic 
current, J, generally decreasing with decreasing p. In the 
partially asymmetric model pc, and at the same time J 
approaches zero in the thermodynamic limit. More pre- 
cisely the stationary velocity v vanishes in a large ring of 
N sites as: 
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and the accumulated distance traveled by the particles, x, 
in time, t, is given by a; '--^ t^^^ , where z is the dynamical 
exponent. 

In this letter we consider the partially asymmetric dis- 
ordered model and will show that a vanishing current 
state can be encountered for pt disorder, too. We will 
derive an analytical expression for Zpt (pt disorder) and 
conjecture a simple relation for Zst (st disorder). We will 
also consider a new stationary state in which the average 
drift is zero and the system has a diffusive motion. We 
show that this new state for both types of disorder has 
an infinite randomness fixed point (IRFP) scenario. So 
far IRFP has only been observed for quantumj^ and 
stochastic models |0, 0| on disordered lattices. 

Here we consider first the ASEP with pt disorder and 
a configuration is characterized in terms of the number 
of empty sites, n^, in front of the ith particle. The sta- 
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tionary weight of a configuration m , n2 
by Plli: In {ni,n2, ■ ■ ■ ,nM 



, um is given 
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provided pi > for all particles. The stationary velocity 
is given by: 



Zn-1 



M 



fN{{n,}) . (3) 



ni,n2,...,nM 



where in the summation = N — M. We are 

interested in the properties of the state in the thermody- 
namic limit, where we define the control-parameter as: 




FIG. 1: Renormalization scheme for particle clusters. If 52 
is the largest hopping rate, in a time-scale, r > 1/52, the 
two-particle cluster moves coherently and the composite par- 
ticle is characterized by the effective hopping rates q and p, 
respectively, see the text. 
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such that for 5 > ((5 < 0) the particles move to the right 
(left). Here and in the following we use [. . . ]av to denote 
averaging over quenched disorder and var[a;] stands for 
the variance of a;. In the following we restrict ourselves 
to the domain 6 > 0. 

The consequences of these formulae for random hop- 
ping rates have been analyzed thoroughly if v is non- 
vanishing in the thermodynamic limit |lfl llll IT^ . in 
particular for the totally asymmetric model. For the par- 
tially asymmetric model we consider the maximal term, 
dmax = ina-xlg^} = which in Eq.Q is dominated by 
a product, Y[a<j<b1j/Pj^ where a < j < 6 is the largest 
region in which qj > pj. The probability of existence of 
this rare region is, P{b — a) ~ exp(— a(6 — a)), therefore 
among M particles its typical size is 6 — a In Af. So we 
obtain: gmax ~ exp(cr(6 — a)) Af, 7 > 0. As a conse- 
quence, in the thermodynamic limit Z^^m is dominated 
by that therm, in which n^. is macroscopic and therefore 
the stationary velocity in Eq. Q is given by: v = g^\^ . 

The distribution of the quantities, g^, for a fixed /i and 
in the thermodynamic limit g^^ = \ is, up to a prefactor, 
in the form of a Kesten- variable [23|. For large A this 
distribution takes the form, P(A) ^ A~^^"'"^/^\ where z 
is the positive root of the equation 



l/z 
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We argue that Zpt = z. For a large, but finite M, and for 
a given realization of disorder two variables, and g^', 
have negligible correlations, provided \fi — fi'\ > ^, where 
^ is proportional to the (finite) correlation length in the 
system. Consequently, the distribution of the 5^ vari- 
ables in a given sample has the same power-law asymp- 
totics with the exponent in Eq.jSJ. Since the stationary 
velocity is the inverse of the largest g^, the distribution 
of V for different samples is obtained from the statistics 
of extremes 21]. Here we use the result that the distribu- 
tion of the maximum of independent random variables, 



which are taken from a distribution with a power-law 
tail is universal and given by the Frechet distribution: 
P(u) = u-i-i/Vze"""'^' with u = gmaxCM-\ Thus 
V ^ Af"^ and with M/N = 0{l) we obtain from Eq. (UJ 
the announced result. For small 5 the dynamical expo- 
nent is divergent: z ~ {25)^^. 

At (5 = 0, when the stationary velocity in the system is 
vanishing and the above formalism does not work, we use 
a strong disorder renormalization group (RG) approach, 
which is analogous to that applied recently for absorbing 
state phase transitions with quenched disorder and 
orig inates in the theory of random quantum spin chains 
|l6l 12^ and random walks the RG method one 

sorts the transition rates in descending order and the 
largest one sets the energy scale, = max({pi}, {g^}), 
which is related to the relevant time-scale, r = fi"^. 
During renormalization the largest hopping rates are suc- 
cessively eliminated, thus the time-scale is increased. In 
a sufficiently large time-scale some cluster of particles 
moves coherently and form composite particles, which 
have new effective transition rates. To illustrate the 
method (see Fig. let us assume that the largest rate 
is associated to a left jump, say ^1 — furthermore 
92 ^ P2,q_i,Pi- In a time-scale, r > the fastest 

jump with rate 52 can not be observed and the two par- 
ticles 1 and 2 form a composite particle. The composite 
particle has a left hopping rate q = qi, since a jump of 
particle 1 is almost immediately followed by a jump of 
particle 2. The transition rate to the right, p, follows 
from the observation that, if the neighboring site to the 
right of particle 2 is empty it spends a small fraction 
of time: r = P2/{P2 + 92) ~ P2/92 on it. A jump of 
particle 1 to the right is possible only this period, thus 
p = pir « PiP2/q2- The renormalization rules can be 
obtained similarly for a large p: 
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The RG scheme outlined above is completely equivalent 
to that of a random antiferromagnetic (dimerized) XX 
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FIG. 2: Scaling of the velocity distribution in the Grif- 
fiths phase at two concentrations (c — 0.2 and c = 0.3) 
for pt (a,c) and st (b,d) disorder. Data denoted by sym- 
bols -|-,x,0,n correspond to M = 64,128,256,512 for (a,c) 
and N = 64, 128, 256, 512 for (b,d), respectively. The Frechet 
distribution with the dynamical exponent: z = Zpt, calcu- 
lated from Eq.J^ and with Zst ~ Zpt/2, as given in Eq.® are 
indicated by a full line. 

spin chain of 2M sites defined by the Hamiltonian: 

2M 

Hxx = -Yl J. {S^St,, + SfSf^,) , (7) 

i=l 

with J2i_i = Pi and J2i — qt- Here 5 in Eq. plays 
the role of the dimerization. The strong disorder RG for 
the random XX-chain is analytically solved |3, and 
the presumably asymptotically exact results|25j can be 
directly applied for the ASEP with pt disorder. In an 
extended part of the off-critical regime, (5 > 0, the cor- 
relation length is finite, but the typical time-scale, tr, is 
divergent. This is the so-called Griffiths phase 05 
which several dynamical quantities, such as the suscep- 
tibility are singular. In the RG procedure in this phase 
almost exclusively the left hopping rates are decimated 
out |2Q]. After Af-steps of decimation we are left with 
a single particle having a vanishing q/p and p ^ Af~^, 
with the dynamical exponent given in Eq. (jSJ . This com- 
pletely coincides with our previous results. At the critical 
point, (5 = 0, left and right hopping rates are decimated 
symmetrically and the system scales into an IRFP. After 
M-steps the remaining effective particle has a symmetric 
hopping probability: q ^ p exp{—constM^^^) . Thus 
the motion of the system is diffusive and ultra-slow, the 
appropriate scaling combination is given by: InvM^^^^ . 
Close to the critical point the correlation length in the 
system, ^, which measures the width of the front, is given 

These analytical results have been checked by calcu- 
lating the velocity distribution of a periodic system with 
M/N = 1/2, using a bimodal distribution with piqi = r, 
for all i, and P{p) = c8{p—l) + {l — c)5{p — r)^ with r > 1 
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FIG. 3; Scaling plot of the velocity distribution at the critical 
point for pt (a) and st (b) disorder. In the latter case we have 
included a logarithmic correction term. 

and 0<c<l/2. In this case the control-parameter is 
(5 = (1 — 2c)/[2c(l — c)lnr] and the dynamical expo- 
nent from Eq. (jSj is z = lnr/ln(c^^ — 1). (In the limit 
c ^ 1/2 one can make a direct calculation by noting that 
an extremely small velocity of v{l) = Cr^' can be found 
in such a sample, in which I consecutive q rates have 
the value of r > 1. Such a sample (rare event) is real- 
ized with an exponentially small probability of P{1) ~ c'. 
Averaging over I we arrive to 2: ~ — Inr/ In c.) As seen 
in Fig. |3ic in the Griffiths phase for two different values 
of the concentration, c < 1/2, the distributions are well 
described by the Frechet statistics and the measured z 
agrees very well with the analytical result. At the criti- 
cal point, c = 1/2, as shown in Fig.|3^, a scaling collapse 
is obtained in terms of the scaling variable, IntiA/^^/^, 
which corresponds to the behavior at an IRFP. 

For st disorder the rate of a jump to i ^ — > i) 
is Pi {qi) and we use the same distribution, as introduced 
before. In the small c limit z^t can be calculated along 
the lines presented for pt disorder, but now the veloc- 
ity in a rare event of length / is [23 v{l) ~ ^r^'/^+^/^. 
The factor 1/2 in the exponential is due to the fact, that 
the boundary between the p = 1 and p = phases lies 
in the middle of the unfavorable domain (particle-hole 
symmetry). For small concentration we obtain: 

Zst = Zpt/2 . (8) 

According to our numerical investigations, as shown in 
Fig. Ebd, the relation in Eq. seems to hold for not 
small values of c, too. In particular, at the critical situa- 
tion, c = 1/2, we recover the result of the IRFP scenario, 
see Fig. |3Jd. 

In the following we present arguments, why the re- 
lation in Eq. ^ can be generally valid. The station- 
ary state both for pt and st disorder has a macroscopic 
phase separation. However, for pt disorder the occupied 
region slowly moves, which is due to single hole diffu- 
sion into the opposite direction, which takes place in a 
position dependent random potential, the hopping rates 
of which, Pi and qi, are generated by the particles. The 
stationary velocity is the same as for a Sinai walker j29l | , 
see Eq. ^ and the relevant time-scale in the problem, 
T ~ , is given by the time needed for a single hole to 
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overcome the largest barrier (rare event) in the sample. 
This occurs with a probability of Ppt(T) ~ r~^/^p*, since 
the typical value of t = Tpt in a large sample is given 
by Ppt{Tpt)N = 1, thus Tpt ^ N^^*, in accordance with 
Eq.©. 

On the other hand for the st disorder the position of the 
occupied block is fixed and the diffusion of holes through 
the occupied phase and diffusion of particles through the 
empty phase will result in the stationary current. In this 
case, due to particle-hole symmetry, the rare event con- 
sists of two (independent) large barriers, one for the holes 
and one for the particles, both having the same time- 
scale. The probability of occurance of this rare event is 
Pst (t) ~ Ppt (r) . Now the typical value of r = Tst is given 
by PstiTst)N = 1, thus Tst ^ N^'p*/^ ^ N^^*, from which 
Eq. ® follows Ha. 

In summary the ASEP with sufficiently strong parti- 
cle or lattice disorder has similar behavior. At (5 = in 
both cases the IRFP scenario holds, whereas for ^ > 0, in 
the Griffiths phase we have a singular dynamical behav- 
ior governed by a dynamical exponent, which is given in 



Eqs. (O and The case of pt disorder can be described 
through RG transformation as an effective single-particle 
problem. For the st disorder case many particle effects 
seem to be important. Our findings might have possible 
applications, e.g. in the case of intracellular transport, 
which typically takes place on one-dimensional tracks. In 
case of uni-directional stochastic motion imperfections of 
the tracks generically lead to local effects, as far as phys- 
iological relevant particle concentrations are considered. 
By contrast we have shown that in case of bi-directional 
motion condensation of particles is generically observed 
in the presence of strong disorder, including the possibil- 
ity of an effective control of the particle velocity according 
to Eq. (HI). 
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